Evolution equation for bidirectional surface waves in a convecting fluid 
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Surface waves in a heated viscous fluid exhibit a long wave oscillatory instability. The nonlin- 
ear evolution of unidirectional waves is known to be described by a modified Korteweg-deVries- 
Kuramoto-Sivashinsky equation. In the present work we eliminate the restriction of unidirectional 
waves and find that the evolution of the wave is governed by a modified Boussinesq system . A per- 
turbed Boussinesq equation of the form ytt - y^x - e^ {yxxtt + {y'^)xx) + {yxxt + yxxxxt + {y'^)xxt) = 
which includes instability and dissipation can be derived from this system. 

^ ■ PACS numbers: 47.20.Bp, 47.35.-i, 47.35.Fg 
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H ' I- INTRODUCTION 

, 

Long wave oscillatory instabilities arise in different contexts. They have been found in the study of unstable drift 
waves in plasmas , fluid flow along an inclined plane 0, Q , convection in fluids with a free surface 0, 0, IE > the 
y—^ Eckhaus instability of traveling waves and, more recently, in solar dynamo waves 0|. The nonlinear evolution of 
the instability is generically described by a Korteweg-deVries-Kuramoto-Sivashinsky (KdV-KS) equation 

I ~t~ I^I^X ~t~ ^^XXX '^'^XX ^XXXX 0; 

in some cases with different or additional nonlinearities B 13 ■ The solutions to this eqiiation have been studied 
' ^ , extensively focusing mainly on the formation and interaction of solitary pulses ITll ll^ . It exhibits parameter 
^ ■ regimes where different behavior is found, numerical studies have shown the appearance of modulated traveling waves, 

■ period doubling and chaos |0, . Exact solutions have been constructed by different methods |3 0, 0, ■ The 
limits of large and small dispersion 5 have been studied rigourously 0, 0| . 

^ The derivation of this equation in the physical situations mentioned above involves the assumption of unidirectional 

■ waves. This is natural in the situations where a preferred direction exists; the electron diamagnetic direction in the 
case of plasma waves, towards the equator in the case of solar dynamo waves, downhill for flow along the inclined 

, plane. However in the case of surface waves in a convectin g flu id this assumption is not natural. The assumption 
of unidirectional waves for surface waves was abandoned in |2Ct l2l| and a modified Boussinesq system of equations 
\^ , derived. 

■ The purpose of this article is to reexamine this problem to find the nonlinear evolution equation of the instability 
\ without the assumption of unidirectional waves, which will be the natural extension of the KdV-KS equation to 

. ^ . bidirectional waves. A higher order Boussinesq system is derived which, following the usual procedure for water waves 
|23l |. can be reduced to a modified Boussinesq equation including instability and dissipation. Our main interest is 
to find the form of this equation, rather than focusing on a particular physical problem. The equation we find has 
• been introduced as one of two possible heuristic extensions of the Boussinesq equation to dissipative systems p^ . 
Jv>( , to our knowledge it has not been derived in a systematic way for any physical problem. In this work we show that 
; I the modified Boussinesq equation is derived systematically for systems with long wave oscillatory instabilities, in a 
_ manner analogous as the Boussinesq equation is derived for surface waves in an ideal fluid. For this reason we take 
as a model the simplest case in which this instability can be found, more precisely, we consider surface waves in a 
convecting fluid neglecting surface tension. 
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II. FORMULATION OF THE PROBLEM 



We consider a layer of fluid which, at rest, lies between 2 = and z = d. Upon it acts a gravitational field g = —gz. 
The fluid is described by the Boussinesq equations 

V -v = Q (1) 

Po^ = -Vp + iN^v + gp (2) 

where the density p depends linearly on the temperature 

p = po[l-a(T-ro)]. (4) 

Here ^ = ^ + w • V is the convective derivative; p, T, and v denote the pressure, temperature and fluid velocity 
respectively. The quantities po and Tq are reference values. The fluid properties, namely its viscosity /i, thermal 
diffusivity k, and coefficient of thermal expansion a are constant. Furthermore we restrict ourselves to two dimensional 
motion so that v — (u, 0, w). 

The fluid is bounded above by a free surface on which the heat flux is fixed and upon it a constant pressure pa is 
exerted. Below, it is bounded by a plane stress-free surface which is maintained at constant temperature. As the fluid 
moves the free surface is deformed, we shall denote its position by z = + 77(2:, t). The boundary conditions on the 
upper surface are 

T?t + U'qx = w (5) 

P - Pa- -J^lWz + UxItI.- 'qxiUz + Wx)] = Q (6) 

^^{'^-vl){uz+Wx) + '2.p.r]^{wz-u,^) = Q (7) 
n • VT = -F/k (8) 

on z = c? + 77. Here, subscripts denote derivatives, = (1 + rj^)^/"^^ h = (— ry^, 0, 1) /N is the unit normal to the free 
surface, F is the prescribed normal heat flux, and k is the thermal conductivity. 

Denoting by Tf, the fixed temperature of the lower surface, the boundary conditions on the lower surface z — Q are 

w = = 0, T = n (9) 

The static solution to these equations is given by = —F{z — d)/k + Tq, ps — po[l + {aF/k){z — d)], and 
Ps — Pa — 9Po[{z — d) + [aF/2k){z — d)^]. We have chosen the reference temperature Tq as the value of the static 
temperature on the upper surface. The temperature on the lower surface is then Tb = Tq + Fd/k. Equations (1 - 9) 
constitute the problem to be solved. We shall adopt d as unit of length, d"^ / k as unit of time, pod'^ as unit of mass, and 
Fd/k as unit of temperature. Then there are three dimensionless parameters involved in the problem, the Prandtl 
number a = p,/poK, the Rayleigh number R = pogaFd'^/kKu, and the Galileo number G = gd^p^/ p?. 

The linear stability theory has been studied elsewhere 0, Q . The result of interest in the present context is the 
existence of an oscillatory instability at vanishing wave number Oc = with critical Rayleigh number Rc = 30. 
The frequency along the marginal curve is given hy u) — aasfG + . . . which vanishes at criticality. In dimensional 
quantities this frequency is that of long surface waves in an ideal fluid. In this problem the damping due to viscosity 
is compensated by heating. 

III. SMALL AMPLITUDE NONLINEAR EXPANSION 

To study the evolution of the oscillatory instability slightly above onset we let 

R = Rc + e'^i?2 

and, since the instability occurs at vanishing wavcnumbcr and the frequency along the marginal curve is of the order 
of the wave number, we introduce scaled time and horizontal coordinates 



^ = ex, T — et. 
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Next we look for a perturbative solution expanding the dependent variables as 

p = psiz) + e^(po + epi + £^^2 + • • •) 
T - Ts{z) + e\6o + eOi + 0^p2 + . . .) 
u = e^(uo + + e^U2 + . . .) 
w — e^(wo + ewi + e'^W2 + ■ ■ ■) 
V = e^(?/o + f^Vi + e^?72 + ■■■), 

and proceed to solve at each order. Our main interest is to obtain a modified Boussincsq system including instability 
and dissipation. This requires that the calculation be carried out to order e^. For the sake of simplicity we choose 
the Prandtl number and the Galileo number equal to 1. This choice, and the neglect of surface tension effects, does 
not introduce any qualitative effects. The details are given in the Appendix, here we quote the relevant results. 
In leading order the horizontal velocity is given by 

where /(^, t) is an arbitrary function to be determined. At each higher order a new arbitrary function appears in the 
solution for the horizontal speed. We will find a coupled system of equations for the surface height and horizontal 
velocity which are a modified Boussinesq system. The first equation of this system is obtained from the kinematic 
boundary condition for the free surface, the second equation from the solvability condition of the equations (1-9). 
The kinematical boundary condition is, at each order, 

VOr +.h=0 (10) 

vir+g^^o (11) 

772. +h^+ af^^^ + (/77o)c = (12) 
V3r + + affcc? - /3/5«r + ifvi + 9Vo)^ = 0. (13) 

The functions h and j are arbitrary functions of ^ and r. The coefficients a and /3 are positive numbers, 
a = 71/168, /3 — 5149/30240. The second equation needed to close the system comes, as mentioned above, from a 
solvability condition. At each order we find, 

VOi +fr=0 (14) 

Vii +9t = (15) 
V2^ + hr- jf^^T + ffi + 30r/o?7oe = (16) 
+ jr - 75«r + rf^^ + p.f^^^^ + i>f^^Tr + (17) 

In these equations r — 2i?2/15 is the measure of instability, the rest are numeric parameters. For our choice ct = G = 1, 
their values are jl = 478/699, D — 4897/10080 and 7 ~ 257/168. Finally we construct a single system, correct to 
order e'^ for the surface elevation r] and arbitrary function F = f + eg + e^h + C'j by simply adding the equations 
above. In doing so some ambiguity arises in the choice of terms of order and e^, where it is possible to replace 
?7j = —Fr + O(e^) or to replace = — ^Tr + O(e^). In the theory of ideal water waves this leads to the existence of 
different families of Boussinesq systems [23. At this stage we use this freedom to minimize the type of terms. With 
this criterion, the kinematic condition, correct to order becomes 

r]r+F^+ ae^F^^^ - Pe^F^^^r + e^{Fr])^ = 0. (18) 

The solvability condition correct to order is 

V^ + Fr- e'7^«r + e\F^^ - t^iji + v)j]^^^r + i^FFi^ + + \(Se^{j]F{)i^ = 0. (19) 

This is a modified Boussinesq system including instability and dissipation in the linear terms of order c" . It describes 
the propagation of bidirectional waves. In the following section, we reduce this system to a single equation for the 
surface elevation, obtaining a modified Boussinesq equation. 
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IV. EVOLUTION EQUATION FOR THE FREE SURFACE 

If we search for unidirectional waves the system (10-17) leads to a modified KS-KdV equation 

Vt + hvVx + S2VXXX + £{rv XX \ ^SVxxxx ~r S4{vVx)x) = (20) 

Our purpose here is to obtain the evolution equation for the surface without this restriction. In this process some 
ambiguity arises as explained in the previous section. In reducing to a single equation we adopt the criterion of 
choosing the form that gives the simplest physically meaningful linear theory. 

A. Linear Theory 

Keeping linear terms in equations H18I19(I . we find that the surface displacement satisfies the equation 

Vtt - Vii - ^^A-q^CTT + e^(»-'7«5T + Bij^i^^r) = 0, 
where A = a + j,ii2=f^ + i^ + l3- Searching for plane wave solutions e'*^e*'^^, we obtain the characteristic polynomial 

5^(1 + e^Ak^) ~ t'{rk^ - Bk'^)s + k^ = Q. 



The growth rate is given by 



and the frequency is given by 



, , le'irk^ -Bk^ 
(s) — 



2 1 + eMfc2 



" =^(^) ^ — tttmp — ■ 

On the marginal curve X — 0, r = ^2k^ . Going back to unsealed variables, we recover the expansion for the Rayleigh 
number R = + a?B valid at small wavenumbers. Small wavenumbers are unstable, large wavenumbers are stable 
as it occurs in the KdV-KS equation. 

B. Nonlinear equation 

Taking cross derivatives of equations (|18I19|I to obtain a wave equation for the surface elevation we obtain 

rirr - ???? - e^Ar^^^rr + e'(r?755. + Sr^^cCCr) + e^[{Fr^\r - \{F^)ii] - 15e'(77')« + 8^Hri^)iir - 0. (21) 

We recognize the order containing F as already appearing in the derivation of the standard Boussinesq equation, 
therefore we only outline the procedure |23l | to obtain the modified Boussinesq equation . We write 

[{Fr^)^-FF^], = -(\r^^ + FA +0{e'), 
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and 

{F\^^2{^l~F^r^)+0{e^). 
Assuming that F vanishes at — cxd we may write 

Replacing these expressions in ()21f) we obtain 

Vtt - Vii - (^Ma-rr + e^(»''7CCr + Brf^^^^r) " 2e^ + J-jri J VrdC'^^ 
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Changing to new variables (Y, T), 

J —oo 

and defining 

N = rj- e^yy^ 

we obtain our main result, a single evolution equation for bidirectional waves including instability and dissipation, 

33 

Ntt - Nyy - (P'ANyytt ~ ye^(^^^)yF + e^(rA^yyT + BNyyyyt) + '&e'{N'^)YYT = 0. (22) 

Up to order we recognize the improved Boussinesq equation. Instability and dissipation are included in the linear 
terms of order e^. The nonlinearity of order is an additional nonlinearity which is characteristic of this problem 
(see Eq. (20)), and which is also found in the study of the Eckhaus instability 

V. SUMMARY 

We have studied the nonlinear development of a long wave oscillatory instability of surface waves. When the 
assumption of unidirectional waves is made the surface evolution is determined by a modified KdV-KS equation. In 
the present study we allow bidirectional propagation of waves and are led to a modified Boussinesq system, where 
instability and dissipation are present at high order. A single equation is obtained for the evolution of the surface 
following the method used in ideal surface waves. In the case of ideal surface waves one is led to a Boussinesq equation 
, here we obtain a modified Boussinesq equation including instability and dissipation. 

Previous studies have dealt with the damped Boussinesq equation (r < 0, B=0) where stable structures, solitons 
and periodic solutions, must be sustained by external forcing [2^ . In the present case r > and it represents the 
destabilizing term. The equation we find has been proposed as one of two possible heuristic generalizations of the 
KdV-KS equation to bidirectional waves. Here we show that it is indeed the adequate bidirectional generalization of 
the KdV-KS equation by deriving it from a physical problem, just as the Boussinesq equation is derived for surface 
waves in an ideal fluid. Numerical and analytical results for this modified Boussinesq equation l|22|) will be the subject 
of future work. 
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APPENDIX: PERTURBATION EXPANSION 

In this appendix the equations to be solved at each order and their solution are given. Terms that vanish have been 
omitted. 

In leading order the equation to be solved are 

UQzz = 

WQz = --"04 

POz = 
&Ozz — ~Wo 

subject to 

^i'o(O) = uo.(0) =0o(O) =0 
UOzil) = 0oz{l) = 
Po(l) = Vo 
VOt = wo(l). 
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The solution is given by 

&o = /j7o(z), pq = 770 
where we have defined To{z) = {z^ — 3z)/6. The kinematical boundary condition is 

VOr + /€ = 0. (A.l) 

At order e the equations to be solved are 

Ulzz = P0€ + UOt 

wu = -ui^ 
Piz = 30 6*0 

Olzz = Oot — Wi 

subject to 

u'i(O) =uu(0) =0i(O) =0 
uu{l) = Ouil) = 
Pi{l)=rii+2woz{l) 
VlT = wi{l). 

The solvability condition uuzdz = implies 

+ /r = 0. (A.2) 

The solution at this order is 

ui=9i^,T), wi=-g^z 

ei=g^To{z) + firTi{z), pi = ?7i + 3O/5P1 (^) - 2/j 

where we have defined Ti{z) = {25z — IQz^ + z^)/12Q and Pi{z) = Tiz{z). The kinematical boundary condition 
becomes 

?7ir+5i = 0. (A.3) 

At order the equations to be solved are 

U2zz = + ^iT - "0€| 

W2z = ~U2^ 

P2z = 30 6*1 + Wxzz - WQr 
02zz — 01t — W2 — ^0{{ 

subject to 

W2(0) =W2^(0) =^?2(0) =0 

P2(l) =V2 + 2wi^{l) + 15r]l - T]oPoz{i) 

V2t = W2{'^) - {uom)i- 

The solvability condition U2zzdz = U2z{^) implies 

m« +9r = 0, (A.4) 

and the solution to this order is 

U2 = t) + h^Uiz), W2 = -hz - fa^Wiz) 

02 = h^To{z) + g^rTliz) + f^rTT2{z) + f^^^T3{z) - zr]of^, 

P2=m + 15r?o + 30ff£Pi(2) + hrP2{z) - 2g^. 
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where we have defined U{z) = z^{39 - IBz"^ + z'^)/24, W{z) = {91z^ - 21z^ + z'^)/168, P^iz) = (25 - 30^V 5^4)/i20, 
T2{z) = {2f>z - IQz^ + z^)/l2Q, Ti{z) = {f,z^ - 180^^ + 1134^^ + 5040^^ - 195750)/6O48O. The kinematical boundary 
condition becomes 

7?2r+/i€ + a/£« + (/??o)€=0. (A.5) 

At order we only need to calculate the solvability condition for ^3 and the kinematical boundary condition. The 
equation for the horizontal velocity is 

subject to 

U3,{0) = 0, U3,{1) = -W^{1). 

The solvability condition yields 

V2i + hr- 7/««T + //« + 30?7or?o4 = 0, (A.6) 

and the solution for the horizontal velocity is 

U3 = t) + gaU{z) + h^rUsiz) + (r?2| + hr + fh + 30r?or/0€)^V2- 

In the expression above [73(2) (-1022^2 + 294z'* - 28^*^ + ^^)/672. 
The kinematical boundary condition at this order 

??3t = ^3(1) + (wo?7i){ + 

yields 

%r + k + ag^ii - Phi^T + {fvi + gvoh = (A-7) 

The pressure satisfies the equation 

P3z = 30 ^2 + W2zz - Wlr + ^0$^ + R2O0 

with 

P3(l) = m + 3077o?7i + 2^2^(1) - ?7oPi^(l). 

Its explicit solution is not indispensable to continue the calculation. Finally, the last equation needed is obtained 
from the solvability condition of the horizontal velocity at order e^. We have 

UAzz = P3C + M3t - M2{{ + (moMi){, 

with 

W4z(0) = 0, ^40(1) = -??oM22z(l) - W2|(l) - mwoizi'^) - 2r?o{(woz(l) - W0{(1))- 
The solvability condition leads to Eq. (17). 
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